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ABSTRACT 


Three-Dimensional linear secondary instability is 
investigated for boundary layers with pressure gradient and 
suction in the presence of finite amplitude Tollmien-Schlichting 
(TS) wave. The focus is on principle parametric resonance 
responsible for strong growth of subharmonics in low disturbance 
environment. Calculations are presented for the effect of 
pressure gradients and suction on controlling the onset and 
amplification of the secondary instability. 

PACS numbers 47.15.cb , 47. 20. Ft 



I . INTRODUCTION 


In the study of boundary-layer control, it is well known that 
suction can completely stabilize a laminar boundary layer • 
Laminarization of the flow by suction, and subsequent viscous 
drag reduction, is the principle and the most effective mean used 
for laminar flow control (LFC) [1]. The overall efficiency of 
an aircraft with LFC is reduced by the power consumed by suction 
systems* Hence, it is necessary to keep the boundary layer 
laminar with the least possible suction. To accomplish this one 
needs to : 1) accurately calculate the stability characteristics 
of the flow, 2) understand the laminar to turbulent transition 
process, 3) evaluate the effect of suction control on some of 
the early stages leading to transition, which is our concern in 
this paper. 

Most of the theoretical work for LFC has relied on linear 
primary stability theory and the e n method, where a limited 
growth of the primary two-dimensional (2D) traveling disturbances 
(TS waves) is accepted such that the power requirements are 
minimized and almost all stages of transition are avoided. The 
effect of suction on the primary TS instability is well 
established and known to be drastic. This work has been 
summerized in the report of Saric [2] which lists most of the 
important papers in the subject. 
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Three dimensionality is known to be a necessary prerequisite 
for transition [26] . A key three-dimensional (3D) phenomenon in 
the early stages of transition is characterized by a strong 
secondary instability of 3D disturbances in the presence of 
finite amplitude two-dimensional (2D) primary disturbances. 
Primary 3D disturbances might be stable or very slowly growing 
in the absence of the 2D waves. This secondary instability has 
been recognized experimentally [3-7] in the uncontrolled 
boundary layer, and observed numerically in the uncontrolled 
[8,9] and controlled [10] boundary layers. Several routes to 
transition in boundary layer have been identified. One route is 
characterized by subharmonic 3D disturbances that are excited in 
the boundary layer at very low level of the TS amplitude. This 
mechanism produces either the resonant wave interaction predicted 
by Craik [11] (C-type), or the secondary instability of Herbert 
[12] (H-type). The other route, the so called K-type (or peak- 
valley splitting), occurs as the TS amplitude exceeds a certain 
threshold value. This represent the path to transition under 
conditions similar to Klebanoff et al.[3]. At sufficiently 
large amplitude one of the routes or a mixture of both will 
appear depending on the disturbance background. 

In this paper, we are concerned with the subharmonic 
secondary instability mechanism as a route to transition that 
seams to be most dangerous in low-disturbance environments, as 
they are found in free flight. We investigate the development 
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of a subharmonic secondary instability in a boundary layer with 
pressure gradients controlled by suction. Our objective is to 
evaluate the effect of suction control on this early stage 
leading to transition. Several questions need to be answered. 
Does suction delay the onset of the secondary instability? How 
sensitive is the growth of the secondary instability to the 
intensity of suction? What is the effect of the initial 
amplitude of the primary wave on this sensitivity? What is the 
effect of slight changes in the pressure gradient? Does the 
effectiveness of suction as a method for delaying transition 
depend on where it is applied, or on its intensity? Finally, if 
the boundary layer would be kept laminar with the least possible 
suction, then should one allow for a limited growth of the 
primary wave, or should one increase suction to fully stabilize 
the primary wave?. 


II. PRIMARY INSTABILITY 

We consider a 2D boundary-layer flow of an incompressible 
fluid with inviscid flow field given by U = U(x) and distributed 
suction given by v = vw(x) at the wall, where x is the streamwise 
direction and y is the vertical direction. The flow is governed 
by the nonsimilar boundary-layer equation, 

A 2 . A 

(l) + ff vi +4( 1- fAj ) -Yfy - 2$( f^ f|«| - fj f^ ) 
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with boundary conditions. 


(2) f*|($,0)-0 , f (3,0) + 2|fj (5,0) - 0 

f«|( $ , V«) = 1 

given in GSrtler variables, 

(3) d| - U € dx d*] - ( U e / fll ) dy 

Then the velocity components u and v in terms of the new variables 
are, 

(4) u - lj.f^ 

A 

(5) v * vw - ( U«/ fzj ) [ f + 2$ f § + *| ( £ - 1 )f^ ] 
and 

(6) Y( l ) - (/2J / U e ) vw(| ) 

(7) £ (1 ) - ( 2$/ U e )(dU e /dx) 


If both the suction and pressure gradient functions are constant 
( equal to and fa respectively), then f( ^ ,*) ) is a function 
of *) only and we have a similar boundary layer governed by. 


(8) 

+ f f rj + fa ( 

1 “ U) ) - 

X‘ri - 0 

(9) 

f(0) = f*j(0) - 0 , 

f»j-* 1 as 

CO 

where 

A 

the condition Y (^ ) 

“ constant 

demands that vw be 

proportional to U e / • 

For the 

case of flat plate vw 

is proportional to 1/ /"x! 
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Stability calculations were carried out for similar suction 
constant and similar pressure gradient = constant using 
the similar equations (8)-(9), and for constant continuous 
suction vw(£) « SL (given by equation (6)) using the nonsimilar 
equations (l)-(2). Negative and negative SL indicate 

suction, while negative indicates unfavorable pressure 
gradient. Equations (8)-(9) were numerically integrated using 
a shooting technique with a Runge-Kutta integrator, while 
equations (1)— (2) were solved using a second-order finite 

difference technique. The mean-flow results were scaled to 
conform with the way stability equations were 
nondimensionalized. 

Now we consider the primary instability of the calculated mean 

flow with respect to 2D quasi-parallel spatially growing TS 

disturbances. Squire's theorem implies that the critical 

disturbance is 2D. Dimensionless quantities are introduced by 

using the reference velocity Ue and the referece length 

»/2 

L * (i*x / U e ) ,so that Reynolds number is given by 
i/z 

R ■ ( U*x /a') , where x measures the distance from the leading 
edge of the plate, and jj is the fluid kinematic viscosity. 

At sufficiently large distance from the leading edge, primary 
instability of the laminar flow occurs with respect to TS 
disturbances. These disturbances take the travelling wave form, 
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(10) (Up.Vp.Pp) = A 0 [ u(y) ,v(y) ,p(y) ] exp [i( Jo< dx - u) t) ] + c.c. 

where for the spatial stability analysis o< is a complex 
wavenumber given by o< =®< r + io<{ and CO is a real disturbance 
frequency, and c.c. denotes complex conjugate terms. The 
eigensolutions u,v,and p are governed by a fourth-order system 
of equations that is given in the Appendix. This system is 
numerically integrated as initial value problem using a 

combination of shooting [ 13 J and Newton-Raphson iteration 
technique that employs a Gram-Schmidt orthonormalization 
procedure. 

The linear stability theory of primary instability provides o< 

for a given w and R. Then the integration of the growth rate 

-off gives the amplification factor (or the amplitude ratio), 

R 

(11) In ( A / A # ) - -2 Jot ; dR 

Rop 

where A 0 is an arbitrary initial amplitude of the primary 
instability at R 0 p (R where the onset of the primary wave) . 

The eigensolutions may be normalized such that A measures 
directly the maximum r.m.s. value of the streamwise disturbance, 
that is 

2 

(12) max | u(y) I -1/2 
o $ y < oo 
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Since the primary instability of boundary -layer flows is induced 
by viscosity, the growth rates and amplification factors here 
are typically very small compared to the convective length 
scale. 


III. SECONDARY SUBHARMONIC INSTABILITY 


The basic state under consideration is the calculated 2D flow 
with suction and pressure gradient at finite amplitude A of the 
primary TS wave, that is 


(i3) ( u w ,v i>, p b> = (u,0,p) +A [ u (y)> v (y)»p(y) ] expue) 


Where 


A 


A 

O 



0 = 


/ 


cv r dx - u) t 


assumed constant, and 


We consider the 3D quasi-parallel spatial subharmonic 
instability of the basic flow given by (13). The finite 

amplitude primary wave acts as a parametric exitation on the 
secondary instability. Following the analysis of Herbert [12,14] 
and Nayfeh [15], we apply Floquet theory and express the 

secondary wave using the normal mode concept, 

(u s ,v s ,p s ) = exp( fVdx)[u(y),v(y),p(y)] exp(~i©) cos/S z + c.c. 
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( 14 ) 


w« 


= exp(J V dx)w(y) exp(-^iQ) sin/5 


z + c.c. 


where ft is a spanwise real wavenumber, and Y = V r + iVj is a 
characteristic exponent. The spatial growth rate of the 

secondary wave is given by Yf , while Yi can be interpreted as a 
shift in the streamwise wavenumber. In our calculations, we 
consider only the case of Yj = 0, that is the secondary wave is 
perfectly synchronized with the basic state. 


The secondary wave (14) is superposed on the basic state (13) 
and the result is substituted into the dimensionless 
Navier-Stokes equations. The meam flow plus the 2D TS 
quantities are substracted, and the resulting equations are 
linearized in the secondary disturbances. Then one obtains an 
eigenvalue problem that can be written as six first-order system 
of ordinary differential equations in the form. 


(15) DZ, - Z t 


(16) DZ Z = Z, + R DU % + R ( Yr+£ ±<* r ) Z g 

+ A R [ ( Y r +j i <* r ) u Z s + v Z z + Du Z 3 ] 


(17) DZ, = - ( Y r +^i »<r) Z 4 - A Z^, 
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(18) 


DZ^ “ Zj 

(19) DZ 5 - ^Z^-r/Zj + A R [ (K r -ji«( r )uy vZ s ] 

(20) DZ g - - R"‘ ( Y*+^i«r) Z 2 - R"V Zj-R^Zs 

+ A { ( Y, -§ i <*V) v Z, - [ ( Y r - i*r) u + Dv ] Z 3 
+ 4 v Z^ } 


where ( ) Indicates a complex conjugate quantity, D»d/dy, 


A * R [ ( V«4i*r) U -4p iU] - ( Yr + xi ) + p , 


A 

Z| “ u 


A 

Du 


A 

Z* ■ v 


A 

W 


Zg * Dw , 


and the boundary conditions are 


( 21 ) 


Z, - Z 3 - Z g - 0 at y* 0 


Z, , Z^ , Z^ 0 as y -*■ £>o 


When A*0,the system of equations ( 15)— ( 21) govern a primary 
subhamonic wave. For Aj*0, this system was numerically integrated 
as initial value problem from y = y e (edge of the boundary 
layer) to the wall. The eigenvalue search used a 
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Newton-Raphson iteration technique to satisfy the last boundary 
condition at the wall. A well tested code SUPORT [13] is used 
which is coupled with an orthonormalization test based on the 
modefied Gram-schmidt procedure to overcome the stiffness of 
the integrated system of equations. For more details on the 
numerical procedure, the reader is refered to El-Hady [16]. 


The linear stability theory of the secondary instability 

provides Y for a given and R. Then the integration of the 
growth rate Yf gives the amplification factor, 

/ 

(22) In ( B / B 0 ) = 2 \ Y r dR 

*oS 

where B o is an arbitrary initial amplitude of the secondary 
instability at R^ (R where the onset of the secondary wave) . 
Since the secondary subharmonic instability originates from a 
strong mechanism of combined tilting and stretching of the 
vortices [18], the growth rates and amplification factors here 
are large and occur on a convective length scale. 


IV. RESULTS AND DISCUSSION 

For the case of no suction and zero pressure gradient, our 
results are in full agreement with those obtained by Herbert et 
al.[14] and by Nayfeh and Ragab [17]. The first authers used 
spectral collecation methods to solve both the primary and 
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secondary stability problems , while the others used a numerical 
technique similar to what is being used in this paper. 

Results are presented in this section to show the effect of 

similar suction parameter , continuous suction SL , and similar 

pressure gradient parameter on the development of the 

subharmonic secondary instability. All results reported here are 

-6 

for the nondimens ional frequency F = <O/r = 60x 10 , that 

remains fixed as a wave of fixed physical frequency travels 
downstream. In the present analysis, we limited our 
calculations to a spesific frequency and to small suction and 
pressure gradient parameters. This was done to satisfy the 

assumptions that are inherent in the approximate theory of 

linear secondary instability, namely the periodicity of the 

basic state and the weak variation of the TS amplitude. We 
note that higher frequences and higher suction rates will 
increase nonparallel effects, violating the periodicity 

assumption of the basic state. While lower frequences and 

higher unfavorable pressure gradients will increase the 
variation of the TS amplitude in the unstable range, violating 
the second assumption. 

-C 

At F*60 K 10 , a primary instability grows between R 0 p = 554 

and R ( p = 1052 (first and second neutral points), reaching a 
maximum amplification factor of A/A c = 41.679 for Blasius flow 
( = 0 , = 0). A broad band of spanwise wavenumbers of 
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primary 3D subharmonic waves are subject to amplification in this 
region, but the time and length scales of this instability is so 
small to bear any resemblance to experimentally observed 
transition. A strong growth of subharmonic disturbance can be 
due to parametrical exitation by the finite amplitude TS wave 
[25]. 


4. a Effect of suction 


At R*1050, growth rates of the secondary instability is 

shown in figures 1 and 2 as function of the spanwise wavenumber 
3 A 

b * 10 p / R, for various amplitudes A of the primary wave. 
Figure 1 shows results for Blasius flow, while figure 2 shows 
results for Y^“-.l, “0. These figures illustrate first of 

all the destabilizing effect of A at fixed F and R. Second, at 
very small amplitudes, considerable growth rates exist in a 
small band of wavenumbers ,that extends to larger values as 
the amplitude increases. Third, the maximum growth shifts 

slightly to occur at higher b as the amplitude A increases. 
Fourth, small suction rates has strong stabilizing effect on the 
subharmonic secondary instability. As R increases, a similar 
increase in the growth rates exist at fixed F and A . 


i 
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At R * 1050, figure 3 shows the effect of various suction 
rates on the secondary growth for A =.01. Increasing 

suction at fixed R and A decreases the growth rate of the 
subharmonic instability and limit the band of dangerous 
spanwise wavenumbers . The wavenumber of maximum growth rate 
is b = .17 and is not affected by suction. 

By specifying the initial amplitude of the primary, say 
A o =.001, we can combine the effect of increasing amplitude A and 
increasing R for various suction rates at fixed F. For 

comparison purposes, the amplification factor of the subharmonic 
is calculated using equation (22) from R eS to any R. All 
results shown here are for spanwise wavenumber b=.15 which is an 
avarage value for wavenumbers of maximum growth at various A . 
Figure 4 shows the variation of growth rates Yf* of the 

secondary wave with R for various suction rates, while figure 5 
shows their amplification factors. The growth rates -°({ and 
amplification factors of the primary wave ( In A/A„) are also 
shown in these figures and indicated by dotted curves. 
Initially, the primary instability sets in at R 0 p and - °<i grows, 
then the secondary subharmonic instability sets in at R oS and 
starts to grow strongly due to the increase in A with 
increasing R. Ultimitly, - <*i decays, while Yr reaches a 

maximum at a location where the amplitude of the primary starts 
to decrease. Small suction rates at fixed A delays the onset 
of the subharmonic secondary instability , and decreases 
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significantly its growth rate as well as its amplification 
factor. Figure 5 shows a reduction of the amplification factor 
of the secondary wave (In B/B 0 ) from nearly 28 to 8 due to an 
increase in from 0 to -.05. Increasing suction rate to -.1 
dampen completely the subharmonic instability, although the 
primary shows some growth. This indicates that the onset of the 
subharmonic instability requires that the primary amplitude 
exceeds a threshold value. Notice that this threshold value 

depends on Reynolds number, it decreases as R increases. For 
example, for * 0 the onset of the subharmonic secondary 
instability is at R“740 with a threshold amplitude A * .0029, 
while for ■ -.05, the onset of the subharmonic secondary 

instability is at R“850 with a threshold amplitude A = .0024. 
For ■ -.1, the primary amplitude reaches a maximum of only 
.0014 which apparantly is below the value needed to induce 
secondary subharmonic instability at R > 850. 

The onset of the subharmonic secondary instability as well as 
the maximum amplification factor it reaches are also dependent 
on the initial amplitude A o . Figure 6 ,a case of suction rate 
lf e =-.05, shows a primary instability that sets in at R ■ 650. 
The onset of the subharmonic instability occurs at R ■ 850 for 

A “.001, at R “ 775 for A “.002, and at R “ 635 for A ■ .0066 
(well before the onset of the primary) . The amplification 
factor reaches a value of 8,14.5, and 30 respectively. When the 
initial amplitude A 0 is large enough, the initial instability 
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can be so strong and secondary instability occurs directly 
by-passing the usual growth of the linear primary instability. 

This phenomenon has been documented in cases of roughness and 
high freestream turbulence. In a situation like this the flow 
quickly becomes turbulent [23] and transition prediction schemes 
based on linear primary theory fail completely. 

Figure 7 shows the effect of suction on the subharmonic 
instability at various initial amplitude of the primary . The 
maximum of In B/B„ is used as a basis for comparison . The 
effect of suction on the primary is also shown in the figure and 
indicated by the dotted curve. The figure suggests that 

secondary subharmonic instability is very sensetive to and can 
be controlled by slight suction rates. 

The onset of the subharmonic secondary instability is one 
important feature of the transition process. For LFC purposes, 
one might try to avoid or delay this instability by using 

suction. Then, one faces the question whether suction should be 
applied before or after the onset of the secondary instability. 
To answer this, we take Blasius case as a basis for comparison 
where R 0 p = ^54, R,p = 1052, R os » 740, R, s = 1250, and maximum In 
B/B 0 ■ 28.5. Continuous suction SL =-.035( ““*05 at R * 1000) 

was applied starting at R * 554 (R 0 p for Blasius flow). As a 
result, the onset of the subharmonic instability is delayed to 
R * 800 and continuous to R =1245. reaching maximum In B/B 0 
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12.5. When continuous suction with the same level SL —.035 


was delayed and applied at R * 740 (R QS for Blasius case), the 
onset of the secondary instability was almost not affected, it 
occurs at R * 750 but the maximum of In B/B w reaches 18.5. 
Doubling the suction level SL —.070 ( Y e — .1 at R - 1000) that 
begins at R * 740 shows no effect on delaying the onset of 
secondary instability but bring the maximum of In B/B 0 down to 
10 . 


Two factors can affect the aforementioned behavior of the 
subharmonic instability ,the initial amplitude A 0 of the primary 
and its evolution. The first factor is held fixed in the 
previous comparison, its effect is given before in figure 6 • 

To explain the effect of the evolution of the primary 

amplitude, we included the amplification factor curves of the 
primary amplitude (dotted curves) in figure 8 together with 

those of the secondary instability for the aforementioned 

cases. The figure indicates that continuous suction resulted in 
delaying the onset of the primary instability allowing the 
primary wave to travel further downstream such that its 
amplitude will reach a threshold value needed to induce a 
secondary subharmonic instability . While in cases where 
suction started at R * 740, the primary instability was not 

delayed but its amplification factor is only enhanced by 
applying suction, resulting in reduced primary and hence 

reduced secondary amplification factors. These calculations 
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show that suction should be applied further upstream near R C p 

to control the development of the primary amplitude and not 

near the onset of the secondary instability. Similar 
conclusions were reached by Reed and Nayfeh [19] and by Saric 
and Reed [20] investigating the effect of suction on primary TS 
waves, that suction should be concentrated not in the region of 
maximum growth but near R . 

4.b Effect of suction and pressure gradients 


At R * 1050, the growth rate Vr of the secondary subharmonic 
instability as function of the spanwise wavenumber b*10 /r for 
various amplitudes A of the primary exhibits the same features 
given before in figures 1 and 2. As an example, figure 9 
shows a case for Y 0 — .1 and .04, that illustrates again the 
destabilizing effect of A at fixed F and R. Notice that 
)^*-.l has a stabilizing effect, while —.04 has a 
destabilizing effect on the secondary instability. 

Figures 10 and 11 show the effect of pressure gradients alone 
on the growth rates and amplification factors of the subharmonic 
instability when combine the effect of increasing A and 
increasing R at fixed F. Also included in the figure (dotted 
curves) the primary instability for comparison. These 
calculations are for A o =.001 and b «.15. Curve b for = -.02 
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reaches a value of maximum In B/B o = 57, while curve c for 

p g = -.04 is estimated to reach a value of 80. The figure 
indicates that small unfavorable pressure gradient is strongly 
destabilizing. Similar results were given by Bertolotti [24] 
studying the secondary instability of Falkner-Skan flows. 
Figures 10 and 11 also show that pressure gradients produce 
small changes in R os but large changes in the maximum of In (B/B e ). 

A similar conclusion was reached by Saric and Nayfeh [21] 
for the primary 2D instability. 

Figures 12 and 13 show the effect of both suction and 
pressure gradients on Vr and In B/B„ at fixed F using A c = .001 
and b =.15. The figure illustrates the sensitivity of the 
secondary subharmonic instability to small suction and pressure 
gradients. Comparing with figure 11, we find that maximum In 
(B/B„ ) goes from 80 for = 0, ^ =-.04 to 9 for = -.1, = 

-.04. 


Figure 14 shows the effect of suction and pressure gradient 
on the secondary subharmonic instability using maximum In B/B e 
as a basis for comparison. The primary initial amplitude for 
these calculations is A= .001 and the spanwise wavenumber 
b = .15. Increasing A c will have a destabilizing effect that 
can be infered from Fig. 7. 


V . CONCLUDING REMARKS 


Previous calculations show that stabilization of the boundary 
layer by active means ( suction) or by passive means ( 
modifying pressure gradients) or by a combination is a very 
sensetive process. Weak suction rates produce strong 
stabilizing effect on the subharmonic secondary instability 
(decrease growth rates , amplification factors and limit the band 
of dangerous spanwise wavenumbers). While weak unfavorable 
pressure gradients produce strong destabilizing effect on the 
subharmonic secondary instability. 

The onset of the subharmonic secondary instability requires 
that the primary amplitude exeeds a threshold value. This value 
is Reynolds number dependent, it decreases as Reynolds number 
increases. 

When the initial amplitude of the primary is large enough, 
the initial instability can be so strong and secondary 
instability may occur directly by-passing the usual growth of 
the linear primary instability. 

For laminar flow control purposes , suction should be applied 
near the onset of the primary instability to control the 
evolution of the primary amplitude and not near the onset of the 
secondary instability. 
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Recent progress in understanding secondary instabilities may 
prompt modifications to transition prediction schemes to rely on 
a secondary instability theory instead of the primary instability 
theory. A modefied e criterion can be reached that envoi ve the 
amplitude of the primary disturbance or a measure of the 
background disturbance. But the means by which freestream 
disturbances enter the boundary layer (receptivity 
problem [22]) will remain the key for any transition prediction 
method. 
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APPENDIX 


The primary 2D eigenvalue problem can be written as, 

DX, = X z 

i 

DXg, - [ i R ( « U - i») ) +« ] X, + R DU Xj+ i« R 
DX 3 = - io< X , 

-> 2 
DX^ = - R i<XX^ - [ i ( « U - ) + R <x ] x 3 

with boundary conditions 

X, = X 3 = 0 at y = 0 

X, , X 3 0 as y -*• <*> 

where 

X v = u , X t - Du , X 3 - v , X^ = p 


Figure 1 Effect of the amplitude A of the primary on the 
subharmonic secondary growth for zero pressure gradient and no 
suction* a) A “.01, b) A “.008, c) A “.005, d) A “.002. 
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Figure 6 Effect of the initial amplitude A 0 of the primary on 
the onset and amplification factors of the subharmonic 

secondary instability at suction rate X =-.05. a) A^ =.001, 
b) A e =.002, c) A 0 =.0066. 
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amplification factors of the subharmonic at various initial 
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Figure 11 Effect of pressure gradient parameter 
primary and subharmonic secondary amplification factors. The 
primary initial amplitude A *.001. a) =0., b) —.02, 
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Figure 13 Effect of suction and pressure gradient parameters on 
the primary and subharmonic secondary amplification factors. 
Same conditions as in Fig. 12 
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